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INTRODUCTION 

In this contribution we will review a new approach to some nonlinear 
dynamical systems which is related to the q-deformation (or other types of 
deformations) of linear classical and quantum systems considered in |l| . The 
main idea of this approach is to replace constants like frequency or mass, etc., 
which are parameters of the linear systems with constants of the motion of the 
system. This procedure produces from the initial linear system a nonlinear 
one and as it was demonstrated in 0, || the q-oscillator of [[|, |§] may be 
considered as physical system with this specific nonlinearity which was called 
q-nonlinearity. At the example of q-oscillator the constant parameter which 
was replaced by the constant of the motion dependent on the amplitude of 
the vibration was the frequency. 

But this approach may be used for many dynamical systems. So, the 
constant masses in Klein-Gordon and Dirac equations or the constant signal 



velocity in the wave equation may be replaced by the dynamical variables but 
these variables are chosen to be the constant of the motion of the dynamical 
system under consideration. Conceptually we address the problem (which 
is not new) if the physical constant parameters like light velocity, Planck 
constant or gravitational constant are constant parameters in reality or they 
may change, for example, in the process of evolution? The same question may 
be put about such characteristics of elementary particle as electric charge or 
mass, say, of electron or fine structure constant which play the role of constant 
parameters in the dynamical equations of the theory. The example of the q- 
nonlinearity of |T[ showed that a linear equation or system of equations may 
be considered as a limit of nonlinear system characterized by a nonlinearity 
parameter which is not important for small intensities but starts to play 
important role for the high field intensities or for large energy densities. 

Taking into account such phenomenon naturally yields the deformation 
of the linear equations of the motion which are transformed due to the de- 
formation into the nonlinear ones but of specific type. The corresponding 
nonlinear system of equations is simply a " reparametrized" initial linear sys- 
tem of equations in which constant parameters are replaced by the constants 
of the motion of the nonlinear system. This specific procedure of nonlin- 
earization of the initial linear system of equations selects a broad enough 
class of integrable dynamical systems which are simple but could describe 
the mentioned above physical phenomena. Also such dynamical systems 
are appropriate to treat the q-deformations and other aspects of quantum 
qroups 0, ||, [0 and quantum qeometry f8[ as influence of different types 
of nonlinearities in the corresponding physical processes. 

Of course, the formulated approach may be extended in the sense that 
one could deform not only initial linear system to make it the nonlinear one 
using method of replacing constant parameters by constants of the motion. 
In principle, one could deform a simple initial nonlinear system containing 
some constant parameters and to transform it into another nonlinear system 
by replacing the constants with the integrals of the motion of the nonlinear 
system. So, it is possible to extend the class of integrable nonlinear systems 
starting from simple integrable nonlinear systems and reparametrizing them, 
i. e. replacing the constant parameters by the constants of the motion. 
Such approach may be usefull, for instance, in q-deforming general relativity 
or Yang-Mills type gauge theories. Up to our knowledge such classes of 
both linear and nonlinear deformed integrable dynamical systems have not 



been studied till now. It should be noted that the method of obtaining the 
integrable nonlinear systems starting from the linear ones is close in spirit to 
the generalized reduction procedure of considering the integrable dynamical 
systems [|], fllO |. 

The goal of our work is to present a general scheme for the described 
procedure of deforming the linear systems with finite number of degrees of 
freedom (next Section). In subsequent sections we also will give some ex- 
amples of the deformation including infinite number of degrees of freedom 
(wave equation). Then the physical consequencies of the q- nonlinearity as 
blue shift effect 0, deformation of Planck distribution formula 0, [|11|], and 
change of a charge form-factor [0] will be discussed. In the last two sections 
deformed Klein-Gordon equation and Maxwell equation in the vacuum will 
be considered. 



SPECIAL NONLINEAR SYSTEMS OBTAINED FROM LINEAR 
ONES 

We consider a carrier space lZ n . With any matrix A we associate the 
vector field Xa = XiAjd/dxj. This association is a Lie algebra isomorphism. 
By using the associative product of matrices, we can associate powers with 
any matrix A A ■ A - ■ ■ A = A k for any integer k. The number of 
independent ones is given by the degree of the minimal polynomial of A. If 
we associate vector fields with any one of these powers we get 

X k = (A k )™x m ^- (1) 

and 

[X k , X,} = 0, (2) 

therefore we get mutually commuting symmetries for our initial dynamical 
system Xa- We notice that for a generic matrix A there will be n- 
independent matrices and they are a basis of all infinitesimal symmetries for 
Xa with coefficients constants of the motion for Xa- 

For nonsingular matrices A k , i. e. det(l + A k ) ^ 0, we can define 
invertible transformations 

^ 1 " A k 



which are symmetries for Xa, i. e. 



Q^AQk = A, 



(4) 



or we can consider e = which are also symmetries for Xa- 

When A is not generic, there are other infinitesimal symmetries which 
are not obtained from powers of A and moreover the algebra of symme- 
tries need not be Abelian. Again we can use either the Cayley map or the 
exponentiation to get finite symmetries out of the infinitesimal ones. 

Now we can make the dynamics and the infinitesimal symmetries nonlin- 
ear by replacing the entries A 1 - e 1Z with constants of the motion for Xa- 
With this procedure, infinitesimal symmetries are no more linear and they 
do not close on a finite dimensional Lie algebra. 

If Ta £ T(1Z n ) is the subring of constants of the motion for X A , and 
Xi, X2, . . . , Xs is a basis of linear infinitesimal symmetries for Xa, we 
get nonlinear infinitesimal symmetries Xf = f x Xi + f 2 X2 + • • • + f s X$ 
for any choice of f 1 , f 2 , ... , f s G Ta- If Xa is " reparametrized" , i. e. 
replaced with / Xa, f € Ta-, previous symmetries are broken except for 
those which preserve /, i. e. l Xj f = 0. 

To make contact with the quantum situation we shall consider previous 
procedure in the framework of Hamiltonian dynamics. If we consider the 
carrier space 7Z 2n , with the symplectic matrix 



the dynamics represented by Xa will be Hamiltonian if t A J+ J A = 0. It 
follows that J A is a symmetric matrix and we can define a function Ja{x) =t 
xJ Ax. This function turns out to be the Hamiltonian function giving rise 
to Xa- In constructing infinitesimal symmetries for Xa by taking powers 
of A, not all powers will give rise to matrices in the symplectic Lie algebra, 
we have to restrict to odd powers, i. e. A 1 , A 3 , . . . , A 2k+1 , . . . , A 2n+1 , 
indeed 




(5) 



\A 2k+l )J + JA 2k+1 = 



(6) 



and we find functions 




(7) 



UB, fc} — f[B,C] 



when =' xBx, we get fj} = 0, therefore odd powers of A give 

rise to constants of the motion which are in involution. For a generic matrix 
A, Xa turns out to be a completely integrable system in the Liouville sense. 

The even powers of A give rise to symmetries which are not canonical. 
We can turn them into nonlinear ones by using the procedure we have already 
illustrated previously. To deal with nonlinear transformations it is better to 
work in a geometrical framework. Therefore on lZ 2n we have the dynamics 
Xa = XiA l jd/dxj, the simplectic structure u = J2idxiAdx i+n and canonical 
symmetries X\, X2, X n associated with odd powers of A and functions 
fi, j'2, • • • , f n , while Y\, Y2, . . . , Y n associated with even powers of 
A and are noncanonical. By exponentiating the noncanonical infinitesimal 
symmetries we get invertible transformations from 1Z 2n into 1Z 2n which 
takes us from one Hamiltonian description for Xa to a different Hamiltonian 
description. 

As a final comment from linear algebra we recall that any matrix A can 
be decomposed into the sum 



where S is semisimple and N is nilpotent, moreover [S, N] — 0. From here 
it follows that the flow associated with A can be written as 



This formula shows that if we want (non trivial) bounded motions associ- 
ated with one Hamiltonian vector field A, we have to impose N = and 
S should have purely imaginary (non zero) eigenvalues. Therefore a generic 
linear Hamiltonian system with bounded motions must be a collection of 
noninteracting harmonic oscillators. To be more specific we turn to consider 
a two-dimensional harmonic oscillator. 



THE HARMONIC OSCILLATOR 

On 1Z A with coordinates ( x, v) we consider the equations of motion 



A = 



S + N 



(8) 



e tA = e ts e tN . 



(9) 




dx k 




k= 1, 2. 



(10) 



We introduce new coordinates 



qk = x k , 



Vk 



k 



1, 2 



and get 



dq k 
dt 

associated with the matrix 



= ^kPk, 



dPk 
dt 



— —^kQk 
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A more convenient basis of symmetries gives rise to the vector fields 

d d 



X 1 



d_ 

d 



Qi ~n Vijr- 



d_ 

d 



Ql 7\ ^~ Pi ~n 

up i oqi 



op 2 dq 2 
9 9 

Y 2 = q2T> \-P2TT- ■ 

op 2 oq 2 



Associated functions are 



u 2,2 u 2,2 



By using Y\ and Y 2 made into nonlinear transformations, we can consider 

x k = fk(h k )Pk, 

Vk = fk(h k )q k . (18) 

This nonlinear transformation is a symmetry for the dynamics but is non- 
canonical, it means 

dx k Ady k = J2 d (fkPk) A d(f k q k ) (19) 
k k 

which is different from J2kdpk A dq k giving rise to standard Hamiltonian 
description of our harmonic oscillator with Hamiltonian 

k 

The picture we are presented with is the following: 
Xa has Hamiltonian description given by (cu, H) or (u), H), i. e. 

ix A u = —dH, ix A & = —dH. (20) 

What happens if we associate a vector field with H by using ul i. e. we 
consider 

i xU = -dH, (21) 

it turns out that X is a reparametrization of Xa by a constant of the 
motion. 

At this point we can consider the coordinate system where u has the stan- 
dard form (Heisenberg coordinates) and if is a deformation of the quadratic 
expression or we can use a coordinate system where H has the standard 
quadratic expression but u is "deformed" (via the nonlinear noncanonical 
transformation). As u gives rise to Poisson brackets ("commutation rela- 
tions") this second choice can be interpreted via a deformation of the com- 
mutation relations and is the starting point for "q-deformed oscillators." 

Introducing complex coordinates in R 4 we consider complex coordinates 



a k = x k + iy k , a * k = x k - iy k . 



(22) 



The complex structure in R 4 is given by the matrix 
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satisfying J 2 = — 1. This matrix defines a complex structure commuting with 
the dynamical evolution 



d 
da k 



— a. 



_d_ 

da", 



)■ 



In the (p, q) coordinates we have new complex coordinates 



and, setting 



n k = a k a k , 
Z,k = fk(n-k)oi k , C = fk(n k )a* k . 



(23) 



(24) 

(25) 
(26) 



Thus transformation is not " analytic" (we notice that analyticity depends on 
the complex structure and here we have two alternative complex structures 
compatible with the dynamics). In these complex coordinates, if we take the 
new point, steming from qantum mechanics, which takes the Hamiltonian as 
a primitive concept for the dynamics, we are naturally led to consider the 
following two Hamiltonians 



in the a-coordinates and 



k 



H 2 = 2 &tk 



(27) 



(28) 



in the ^-coordinates. To compare, we express them in the same variables to 
find 

H^Wn k (29) 



and 

H2 = \Y,fl^k)n k . (30) 

k 

We now use the same bracket for them, say 

{a k} a* k } = i. (31) 

We obtain two different dynamical systems, where the one associated with 
H2 is even not necessarily isotropic. The evolution goes from periodic orbit 
to orbit whose closure is a 2-dimensional torus, i.e. the associated systems 
are completely different. Of course, there is no contradiction. Indeed to have 
the same dynamics we should use different Poisson Brackets, namely, the one 
we get transforming Eq. (31). It is now clear that under quantization these 
complex coordinates go into creation and annihilation operators. Therefore 
for the corresponding commutators we can repeat what we have said for the 
Poisson Bracket. 

As a final remark, we notice that the notion of coherent states is a kine- 
matical notion. The coherent states are discussed in quantum mechanics 
from the point of view of the states the closest to classical one. Accord- 
ing to this physical property the coherent states of the electromagnetic field 



harmonic oscillator were introduced in [ 13| . The coherent states (classical 
Gaussian packets) of a charge moving in magnetic field were introduced in 
l4fl , [fl5f1 . Coherent states of spin were introduced in JIB] . 



We start with a Poisson bracket compatible with a complex structure, we 
consider a canonical chart to be anyone such that 

{Xk,Xj} = iSkj, (32) 

then we associate creation and annihilation operators with these variables 
and define the associated coherent states. It is clear therefore that the ingre- 
dients needed for coherent states are: 

1) a symplectic structure 

2) a compatible complex structure 

and together they imply the existence of an Hermitian structure. When 
there are various structures compatible with a given dynamical evolution, 
we find various coherent states (i.e. associated with various symplectic struc- 
tures) preserved by this evolution. The discussed complex structure has been 



used to construct Jordan-Schwinger map for some functional groups in |T7 



CLASSICAL q-OSCILLATOR AS " REPARAMETRIZ ATION 



In previous section we described the general scheme of deformations. Now 
we give in detail the example of q-oscillator. Also, the q-oscillator math- 
ematical properties and its modifications have been used to associate with 
some physical phenomena [|18|-||26"1. The q-oscillators are interesting from 



the point of view of mathematical structure PTfl-|[32| . Let us start from the 



considering the one-dimensional linear harmonic oscillator with frequency uo 
and corresponding nonlinear q-oscillator following the approach of Ref . |l| . 
The linear second order equation for the coordinate of this oscillator x is 



x + u x = 0. 



(33) 



In the form of the first order equations for the position and momentum this 
harmonic oscillator is described by the system 



P 

x = p 



—uj 2 x, 



(34) 



The variables x and p are real ones. If we introduce complex variables 



a 



1 /- /— i \ 
-=(y/wx + -=p), 

V2 yJU) 



a 



\/2 



P) 



(35) 



the harmonic oscillator is described by the first order differential equations 
system 



a 

• * 

a 



luja . 



(36) 



Both systems (34) and (35) may be rewritten in matrix form of Ref. []J 

X = AX. (37) 



The real two-vector X_ e R 2 and the matrix A for the system (34) are 

The complex vector X and matrix A for the system (35) are 

According to our general procedure let us replace the constant matrix ele- 
ments A\ of matrix A by functions of the constants of the motion. The 
system obtained is nonlinear and it may be integrated by exponentiation as 
easily as the initial linear system. We illustrate this procedure again on the 
harmonic oscillator equations of motion (39) since in this case the matrix A 
is diagonal one. Thus the constant frequency ui should be replaced by the 
constant of the motion u)(a, a*) and we have new system for one and the 
same variable a 

a = —icu(a,a*)a, 

a* = i£u(a,a*)a*. (40) 

If one takes the constant of the motion u(a, a*) to be integral of motion 
independent explicitly on time we have the condition for this function 

_ W(Q ., Q . ) = _« + _« =0 (41) 

which after using the equation of motion (40) may be integrated and the 
solution for this equation is 

u(a,a*) = Sl(aa*), (42) 

i. e. any function Q of the modulus of the amplitude is the frequency which 
is the constant of the motion. We will take this function in the form 

n = cof q (aa*). (43) 

The function f q (aa*) will specify the q-deformation of the harmonic oscil- 
lator and we will take it in the form 

f q (z) = — cosh Xz, q = e\ (44) 



where A is the parameter of q-nonlinearity of the vibrations. In case A — > 
the function f q {z) — > 1. The relation to classical q-oscillator (and its 
quantum counter part) consists of the existence the nonlinear noncanonical 
transform in the phase space of the oscillator which in the variables a looks 
las 

/ sinh Aaa* „ / sinh Acta* „ 

V aa* sinh A ' 9 V aa* sinh A 

The Poisson Brackets of the variables a q reproduce the structure of the 
annihilation and creation operators and the commutation relations of the 
q-oscillator |5) (see, below). The dynamics of the physical variables a due 
to the Hamiltonian 

H = ua*a q , (46) 

which is taken in form to be the same as the Hamiltonian of the harmonic 
oscillator 

H= ua*a } (47) 

is the nonlinear vibration corresponding to Eq. (40). The linear dynamics of 
the variable a corresponding to the Hamiltonian (47) is given by the solution 
to the Eq. (36). It is nessesary to point out that in the suggested approach 
the physical meaning of the variables a as related to physical position x 
and physical momentum p in the sense of the measurement procedure by 
the formula (35) is preserved for nonlinearized dynamics. It means that the 
dynamics of the coordinate x and momentum p of the nonlinear q-oscillator 
is described by the system of equations 



x 



f q (aa*)p, 



p = -uj 2 f q (aa*)x (48) 

where 

aa * = -( ux 2 + I p 2 ). (49) 
2 ijj 

For the deformed equations of motion of q-oscillator as the formula (48) 
shows the momentum is the function of the velocity and position. This 
function may be obtained as the solution to the functional equation 

sinh A x 

p(x, x) = — —j- — — T - T — (50) 

A cosh[^{c<JX / + ± p z {x, x)\\ 



considered as the implicit formula for the giving momentum as the function 
of the position x and velocity x. The solution to q-oscillator equation of 
motion 

a = —iua — — - cosh \aa* (51) 
smh A 

is 

a(t) = ao exp[— iut — — - cosh Acto^o! (52) 
smh A 

where 

ao = a(t = 0) 

is the initial complex amplitude of the nonlinear q-oscillator. The solution 
to the equation of motion for the coordinate x of the nonlinear q-oscillator 

* + " 2 — ^rr cosh 2 {-^[xV + p 2 (x, x)}} = (53) 
sinh A 2uj 

where the function p(x, x) is given explicitly by the relation (50), may be 
written in the form 

/ \ r r ^^^t i r ^ r 2 2 2/ * \ 1 1 1 

x(t) = —{exp[^-^cosh{ — [x iu+p{xo,x )\}\ 
+ 6Xp ^^nm^ COSh ^^ 2 +p2 ^ ' X °^^ 

Xq sinh A i -1 r ^ r 2 2 2/ • \n 

+ 2i\iu ° ^2u? X ° U +P X °'" 

X ^ 6XP ^sTnnT COSh ^2^^ 2 +p2 ( x °' io ^] 

~ e M-^^cosh{^-[x 2 u 2 +p 2 (x , x )]}]}. (54) 

Here xq = x(t = 0) and x = x(t = 0) are the initial position and velocity 
of the nonlinear q-oscillator. In the limit A — > we have the usual solution 
for the linear harmonic oscillator. 

One can find for small nonlinearity A C 1 the approximate expression 
for the momentum solving the equation (50) by iteration method. We have 

P = x[l + ^-^(ux 2 + ^) 2 ]. (55) 



This formula may be interpreated as the negative shift of the mass of the 
oscillator by the factor depending quadratically on the energy of the oscilla- 
tions. 



TWO-DIMENSIONAL CLASSICAL q-OSCILLATOR 

The dynamics of generically deformed two-dimensional oscillator is de- 
scribed in previous sections. In this Section we will give the formula for 
specific q-deformation of this oscillator If one considers two degrees of 
freedom, the amplitudes of the first harmonic oscillator a + and of the sec- 
ond oscillator a_ may be q-deformed by two different methods. One is to 
have frequency of the first oscillator to be dependent only on the energy of 
this oscillator n + = \a + \ 2 and the frequency of the second oscillator to be 
dependent only on its energy n_ = |ct_| 2 . Another method of deformation is 
to make the frequencies of both oscillators to be dependent on the full energy 
of vibrations 

n = n + + n_. 

We now will describe the q-deformed variables a q± in the case of q-deformation 
when the frequency of vibrations depends on full energy. Then we have the 
following non-zero Poisson Brackets 

r * i • * (An) coshnA — sinhnA 

{a q± , a qT } = ta±a T n 2 sinhA , (56) 



and 

Here 
where 



% n^ 

{a„±, a*j \ = r~r[(l ) sinhnA + An± cosh nAl. (57) 

q n sinh A n 

a q± = a±F{n) (58) 



„. . / sinh Xn , . 

f <") = V^Ta- (59) 

The introduced deformation of two harmonic oscillators implies the interac- 
tion of these oscillators which exists due to q-nonlinearity. Thus we have 
not only selfinteraction of the different modes but the mutual influence of 
the motion of one oscillator onto the other. 



DEFORMED WAVE EQUATION 



Now we consider the case of the system with infinite number of degrees 
of freedom. An example of such a system is the wave equation with the 
constant parameter which is wave velocity. This constant parameter we take 
to be unity. Thus we start from the wave equation of the form 

= 0. (60) 

In order to clarify the deformation procedure we represent this equation as 
system of equations for decoupled oscillators. To do this we rewrite this 
equation in momentum representation 

<p(k,t) + k 2 (p(k,t) =0 (61) 

where the complex Fourier amplitude 

(f(k,t) — — ( tp(x, t) exp(— ikx) dx, (62) 

2,71 J 

plays the role of new coordinate. Since ip(x, t) = (p*(x, t) we have 
<p(k, t) = Lp*(—k, t). Eq. (61) describes a two-dimensional oscillator with 
equal frequencies for both modes labelled by k and — k. Writing the Eq. 
(61) in the form 

<p(k, t) = n(k, t), 

n(k, t) = -k 2 p{k, t) (63) 

we have the equations of the form (38) in which frequency u 2 = k 2 , p — > 
7i(k, t) and x —> tp(k, t). Due to this we can deform this linear system 
taking the integral of the motion 

n = /^^ 2 M 2 (M) + |i^l 2 ]} (64) 

in which the function playing the role of complex momentum of A;— field 
mode is solution to the infinite system of equations 

VKM) = F k f q { J dk'^-[k' 2 \<p\ 2 (k',t) + \F k ,\ 2 }}. (65) 



The function f q is given by the Eq. (44). Then the parameter fj, plays the 
role of the initial number of vibrations corresponding to given Cauchy initial 
conditions. 

Thus we made the deformation by method used in the previous Section 
for two-dimensional oscillator introducing the interaction of modes through 
the parameter /i. The deformed wave equation may be written as 

<p(x, t) = -J k 2 f;^Mk, ty kx dk. (66) 

This equation may be rewritten in the form for which only scalar field in time- 
space coordinate (p(x,t) is present. For this we replace in the integrand the 
Fourier components ip(k, t) by its expression in terms of ip(x, t) as well as in 
the integral of motion /j, which, in principle, may depend on the momentum 
vector k through the dependence on the Fourier components of the initial 
field <p(x, 0) and <p(x,0). Thus we have the deformed wave equation 

<p(x,t) — —— J J k 2 exp[ik(x — x')]ip(x', t)f 2 {/i[ip(x, 0), ip(x, 0)]} dk dx', 

(67) 

We have pointed out here that the integral of motion \x is the functional of 
the field ip(x, t) of the special form. Being the integral of motion it depends 
only on initial values of field and field velocities. 

Since the parameter \x is chosen as the common integral of motion for 
all field oscillators it does not depend on wave vector k, and the function 
fq(li) can be considered as a commom factor. Then the wave equation may 
be rewritten in the form 

<p(x,t) = fZ(jM)-j]^<p(x,t). (68) 

We have the differential-functional equation which looks like usual wave 
equation with the wave velocity f q {fJ.) which is constant of the motion. Thus 
the procedure of deformation yields us the nonlinear equation for which the 
velocity of wave propagation depends on the initial configuration of the field 
and its time derivative. This equation may be appropriate to describe a field 
behaviour in a media with strong nonlinear response of its properties to the 
presence of the quanta of the field. 

Now we will obtain the solutions to nonlinear deformed wave equation 
(68). The parameter ji behaves as constant for any choice of the initial 



conditions <p{x, t = 0), (p(x, t = 0). Due to this the solution to the nonlinear 
mode-vibration equation (61) is 



ip(k,t) = ^(p(k,Q){exp[i\k\f q {^)]t + esxp[-i\k\f q {fi)]t} 

+ iyfo 0){exp[z|A:l/ g (//)]t - exp[-i|A:l/«fr)]t} J ^ , (69) 

and 

0) = — / y?(x, 0) exp(— i/cx) dx, 

27T J 

(p(k,0) = — / 0(x, 0) exp(— ikx) dx. (70) 

It is multimode generalization of the solution (54) of the one-mode nonlinear 
oscillator. Thus given initial condition (p(x, 0), <p(x, 0) implies that (p(k, 0) 
and ip(k, 0), and /i are given, too. In terms of these values we have the 
A;-th mode solution ip(k,t) and the solution to nonlinear q-deformed wave 
equation (68) are given by Eqs. (62), (69). It is easy to prove that the 
q-deformed wave equation (68) has the soliton-like solutions 

<p ± (x, t)= <f>(x±f q (ji)t) (71) 

where $ is an arbitrary function. In fact, discussed q-deformation implies 
the existence of nonlinear interaction among the modes. The generalization 
to the case of three space coordinates may be performed following the same 
reparametrization procedure. The q-deformed Klein-Gordon equation is con- 
sidered by this method in |53]]. From the point of view of deformed Poincare 



symmetry group the relativistic equations are discussed in 134] 



QUANTUM q-OSCILLATOR 

To make more clear the relation of the suggested approach to classical 
equations of motion with standard quantum q-oscillator formalism of we 
review the description of the oscillator given in PJ. Let us introduce the 
usual creation and annihilation oscillator operators a and a) obeying bosonic 
commutation relations 

[a,a f ] = l. (72) 



Below we assume the classical dynamical variables to which a and a) cor- 
respond to oscillate with a frequency uj = 1. It is known that the operators 
a, a 1 ", 1 form the Lie algebra of Heisenberg-Weyl group. So, the linear har- 
monic oscillator may be connected with the generators of pure Heisenberg- 
Weyl Lie group. In view of the commutation relation (72) the usual scheme 
for generating the states of the harmonic oscillator is based on the properties 
of the Hermitean number operator n 



a 



[a,n] = a, [a ,\h] = —a. 
Thus constructing the vacuum state |0) obeying the equation 

a|0) = 0, 

and the excited states 



,t n 



\n) 



|0> 



which are eigenstates of the number operator h 

h\n) = n\n),n G Z + 



(73) 



(74) 



(75) 



(76) 



the matrix representation of the operators a and a* in the basis (75) have 
the known expressions 






Vi 
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V2 








V 



(77) 



while the number operator n is described by the matrix 

/ ... ^ 
. _ 10...'' 
n ~ 2 ... 

V / 



(78) 



The Hamiltonian for such a system is defined as 



H 



a^a + aa 



(79) 



The q-oscillators may be introduced by generalizing the matrices (77) and 
(78) with the help of the q-integer numbers n q , 



n q = 



sinh nX 
sinh A 



q = e 



(80) 



Here A and q are dimensionless c-numbers, which appear at this purely 
mathematical level. When A = 0, q — 1 and the g-integer n q coincides 
with n. Then, replacing the integers in (77) and (78) by g-integers we 
obtain matrices which define the annihilation and creation operators of the 
quantum q-oscillator, 



a„ = 



a' = 



( ^ 
^/T q 


v 

/ 
VT q 
^T q 

v 



/ 








2 







V 



since the action of n q on eigenstates \n > is given by 



n„\n >- 



sinh nX , 



\n > 



sinh A 

The above matrices obey the commutation relation 



(81) 



(82) 



[a q ,n] 



"q ■ 



[aj,n] 



,t 



(83) 



but the commutation relations of the operators a q and a q do not coincide 
with the boson commutation relations. Eq. (72) is replaced by 

K,o}] = F(n) (84) 

where the function F(n) has the form 

_ sinh A(n + 1) - sinh Xn 
Hn) ~ sinh A • 

For A = (84) reduces to (72). In addition to the above commutation 
relation there exists the reordering relation 

a q a\ - qa\a q = q~ n (86) 

which usually is taken as the definition of q-oscillators. 

It is worthy noting that the operators a q and a q can be expressed in 
terms of the operators a and a, (see, for example, [|TJ) 



> I q - - MJ \ll I. ' ' ' 

where 



a q = af(h), a\ = f(n)a ] (87) 



m = \/-- 

V n 

The comparison of formulae (45) and (87) shows the complete analogy of 
the quantum q-oscillator and the classical q-oscillator discussed in previous 
sections. We have also 

h q = a\a q (89) 

and 

[a q ,n q ] = F(n)a q , [ a l,n q } = -a\F(h). (90) 

In the Schrodinger representation the evolution operator of the harmonic 
oscillator 

[at a + aat)t 



U(t) = exp 



-no 



(91) 



gives the possibility to find out explicitly linear integrals of motion which 
depend on time 

A{t) = U(t)aU-\t) = e iuit a, 

A\t) = U(t)a^U- 1 (t) = e- lU]t a. (92) 



The matrices of the integrals of motion (92) in Fock basis may be obtained 
from the equations 



A{t)\n) 
A\t)\n) 



Suit 



Vn\n - 1), 



e- iujt VnTT\n + 1). 
Let us now introduce the Hamiltonian 



H = u 



for which the evolution operator takes the form 



U q (t) = exp 



— iut 



(lq(lq -\- QqQjq} 



We have for the integrals of motion 

AS) = U q {t)a q U-\t), A\{t) = U q {t)a\U-\t) 
the following explicit matrix expressions p| 



A q (t) 



Al(t) 



y/Tqe^ 1 "-^ 








v^e 






-i{2 q -l q )uit 



\ 





l(2q — l q )ult 

...x 
... 
... 

/ 



(93) 



(94) 



(95) 



(96) 



(97) 



These operators are the generalizations of the linear integrals of motion (92) 
to the case of nonlinear Hamiltonian. This result is a generalization for 
q-oscillators of such integrals of motion for usual and parametric quantum 
oscillator which have been found in |36|] and discussed for constructing co- 
herent states in [37| and 



ANALOGY OF CLASSICAL AND QUANTUM DEFORMATIONS 



By using the example of the oscillator we clarify now the analogy which 
exists in the suggested approach to deformation of classical systems and 
quantum ones. Now we recall what has been done in the previous Section 
for q-deformed quantum oscillator to make clear the connection with the 
procedure of deforming the classical oscillator of first sections. 

Equations of motion for the harmonic oscillator amplitude a we rewrite 
for another operator 

A = h(a\a)a 

where h is a real function, and its hermite conjugate A] = at h{a]a) in 
the same form 

A = -iuA, A\ = iujA\. 
We have in our Hilbert space the vacuum state \l/o which satisfies 

a^ = 0, A^ = 0. 

Thus we could construct two bases in the vector space. One is the standard 
basis 

at" 
vn! 



which is orthonormal in standard scalar product 

<^ n \^ m >= 6 nm . 
Another basis is constructed using the operator A] 



We define new scalar product in the same vector space which is given by 

< ^ n |^ m >= 5 nm . 

The adjoint with respect to this new scalar product need not coincide with 
the old one. We can define the operators 



b*^ n = Vn + ltf n+ i, b^ n = Vn^ 



n-l 



where * means the adjoint in the new scalar product. These operators satisfy 
the commutation relations [b, b*] = l. Taking the Hamiltonian H = uob*b we 



have for the operators b, b* the equation of motion of the harmonic oscillator. 
Thus for one and the same vector space we have possibility to introduce two 
Hilbert space structures. As for the dynamics we have, like in the classical 
case, two different descriptions. Very much as we did for the classical case 
we can use the new Hamiltonian and the old commutator relations to get a 
" deformed" dynamics. As for the partition function, similarly to the classical 
case, we can use the trace defined via the two different scalar products to get 
the same result, i. e. the partition function depends only on the dynamics 
and not on the particular Hamiltonian description we use. Either we change 
the Hamiltonian and for the old scalar product we obtain new dynamics. Or 
changing Hamiltonian and simultaneously the scalar product we obtain the 
same dynamics. For partition function in such case we have the same value 
that was for nondeformed oscillator. 



DEFORMED PLANCK DISTRIBUTION 

In this Section we will discuss what physical consequences may be found if 
the considered q-nonlinearity influences the vibrations of the real field mode 
oscillators like, for example, electormagnetic fields ones or the oscillations of 
the nuclei in polyatomic molecules. First of all this nonlinearity changes the 
specific heat behaviour. To show this we have to find the partition function 
for a single q-oscillator corresponding to the Hamiltonian H = n q 

oo 

Z(T) = $>xp(-/?n g ) (98) 

n=0 

where the variable (3 is the function of the temperature T~ l . The evaluation 
of the quantum partition function of the q-oscillator yields for the specific 
heat that it decreases for T — * oo as 

Thus the behaviour of the specific heat of the q-oscillator found in |lj] for 
A 1 is different from the behaviour of the usual oscillator in the high 
temperature limit. This property may serve for an experimental check of the 
existence of vibrational nonlinearity of the q-oscillator fields. 



q-Deformed Bose distribution can be obtained by the same method start- 
ing from the Hamiltonian H = §{at, a q }+ and one obtains fl| 



- ^ 
<n>= n -(j— 

6 



\(^) - (n)g) + 2((n3) - no(^)o) + (n*) - no(n s ) ) 



(100) 

in which n is the usual Bose distribution function and 

(n fc ) = 2 sinh ^ £ n fe e -/3(n+i/2) _ ( 1Q1) 

2 n=0 

Calculating the partition function for small q-nonlinearity parameter we 
have also the following q-deformed Planck distribution formula 

< n >= r—r-= A — ; z 777* T-a • (102) 

It means that q-nonlinearity deforms the black body radiation formula 0. 
One can write down the high and low temperature approximations for 



the deformed Planck distribution formula [11]. For small temperature the 



behaviour of the deformed Planck distribution differs from the usual one 

< n > -no = -A 2 ^e-^/ feT . (103) 

For the high temperature the nonlinear correction to the usual Planck dis- 
tribution also depends on temperature 

< n > -n = -6A 2 (^)" 3 . (104) 

As it was seen, the discussed q-nonlinearity produces a correction to Planck 
distribution formula and also this may be subjected to an experimental test. 

As it was suggested in the q-nonlinearity of the field vibrations pro- 
duces blue shift effect which is the effect of the frequency increase with the 
field intensity. For small nonlinearity parameter A and for large quantity of 
photons n in a given mode the relative shift of the light frequency is 

5w_A% _1 



This phenomenon of possible existence of the q-nonlinearity may be essential 
for the models of the early stage of the Universe. 

Another possible phenomenon related to the q-nonlinearity was consid- 
ered in [12[ where it was shown that if one deforms the electrostatics equation 



using the method of deformed creation and annihilation operators the form- 
factor of a point charge appears due to q-nonlinearity. 



NONLINEAR KLEIN-GORDON EQUATION 

To demonstrate how the q-nonlinearity may appear in Klein-Gordon 
equation we start from the consideration of usual Klein-Gordon equation 
with mass equal to zero ( c = 1) 

(— -AMx,i) = 0. (105) 

Let us take the plane wave solutions of the equation, i.e., we represent the 
field (p(x.,t) in the form 

<p{x.,t) = J <p(k, t) exp(ikx) dk, (106) 

where Fourier amplitude 

f(Kt) = J </?(x,t)exp(-ikx) dx, (107) 

plays the role of new coordinate. It satisfies the integral equation 

J </5(k, t) exp(zkx) dk = J (-k 2 )ip(k, t) exp(ikx) dk. (108) 

This integral equation is equivalent to the differential one 

</5(k, t) + k 2 ip{k, t) = 0, (109) 

which is an infinite system of decoupled oscillators with frequencies u 2 = k 2 . 
According to suggested procedure we replace this equation by 

<p(k J t) + k 2 f 2 ( f iMk,t) = o. (no) 



Here the new frequency of k-th mode appeared 

co 2 = k 2 f^), (111) 

where the choice of parameter \x determines the deformation. According 
to our ideology it is possible to take it to be an integral of motion of the 
Klein-Gordon equation. There exists a common integral of motion which is 
the full number of the scalar field quanta, and we take it to be equal ji. As 
well as for the wave equation case the parameter \i behaves as constant for 
any choice of the initial conditions <£>(x, t — 0), <£>(x, t — 0). Due to this the 
solution to the nonlinear mode-vibration equation is the following 

<p(k,t) = ^(k,0){exp[ i \k\f q (fi)}t + exp[- i \k\f q (fi)]t} 

+ ly(k, 0){e W [i\k\f q (f,)]t - exp[-i\k\ f q (»)]t} im 

Here 

1 r°° 
1 r°° 

0(k,O) = ^-py_ oo ^(x,0)exp(-2kx) tbc (113) 

Thus the given initial conditions y?(x, 0), 0(x, 0) determine also </?(k, 0), 
(p(k, 0), and /z. In terms of these values we have the k-th mode solution 
(p(k,t) and the solution to nonlinear q-deformed Klein-Gordon equation 

(|-/;WAMx,t) = 0. (114) 
The constant of motion plays the role of signal velocity. 



q-DEFORMED ELECTRODYNAMICS 

We will consider the system of usual linear Maxwell equations in vac- 
uum for the field E(x, t), H(x, t) expressed in terms of the potentials 
A(x,t), (f(x,t) 

1 d A dip 
c dt <9x' 

H = rot A, (115) 



in the form ( c = 1) 

(^-AMx,t) = 0, (116) 

(J^-A)A(x,i) = 0, 

^ + divA = 0. (117) 
ot 

The suggested approach for scalar field of previous Section may be applied 
if one reexpresses the fields </?(x, t), A(x, t) in Fourier basis 

</?(x, t) = j </?(k, t) exp(ikx) dk, 

A(x,t) = y A(k,t)exp(ikx) dk. (118) 

Then we have 

<f(k,t) + k 2 ip(k,t) = 0, 

A(k,t) + A; 2 A(k,t) = (119) 

together with gauge constraint 

(p(k, t) + ikA(k, t) = 0. (120) 

The approach is based on the taking into account that the equations describ- 
ing the linear forced oscillators of the electromagnetic field are deformed due 
to dependence of the frequency of the oscillations on the energy of the elec- 
tromagnetic vibrations. As in the case of the scalar field we will introduce 
the total "number of quanta" of the defiormed electromagnetic field fi which 
is integral of motion. Then in this case we have the deformed nonlinear 
equations of vibrations 

Cp(k,t) + f q (n)k 2 V (k,t) = o, 

A(k,t) + / g 2 (/i)A; 2 A(k,t) = 0, (121) 

The sense of the function f q (n) is the signal velocity.Due to this we have to 
reparametrize the gauge condition to become 



f q - 1 (nMk,t)+ikA(k,t) = 0. 



(122) 



After this we could reconstruct Maxwell equations in space—time analogously 
to the wave equation case. It should be noted that parameter \x in the Klein- 
Gordon and Maxwell equations is given by Eq. (64). The limit of electrostat- 
ics in the suggested type of q-deformation coincides with the usual electro- 
statics described by the Laplace equation (compare with |]12j] ). It would be 
interesting to take into account interaction of the deformed electromagnetic 
field with the sources, but it will be taken up elsewhere. 
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